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Abstract. In this paper, we address the stability of transport systems and 
wave propagation on networks with time-varying parameters. We do so by 
reformulating these systems as non-autonomous difference equations and by 
providing a suitable representation of their solutions in terms of their initial 
conditions and some time-dependent matrix coefficients. This enables us to 
characterize the asymptotic behavior of solutions in terms of such coefficients. 
In the case of difference equations with arbitrary switching, we obtain a delay- 
independent generalization of the well-known criterion for autonomous systems 
due to Hale and Silkowski. As a consequence, we show that exponential sta¬ 
bility of transport systems and wave propagation on networks is robust with 
respect to variations of the lengths of the edges of the network preserving their 
rational dependence structure. This leads to our main result: the wave equa¬ 
tion on a network with arbitrarily switching damping at external vertices is 
exponentially stable if and only if the network is a tree and the damping is 
bounded away from zero at all external vertices but at most one. 
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1. Introduction. Dynamics on networks has generated in the past decades an 
intense research activity within the PDE control community [2,7,11,14,11’]. In par¬ 
ticular, stability and stabilization of transport and wave propagation on networks 
raise challenging questions on the relationships between the asymptotic-in-time be¬ 
havior of solutions on the one hand and, on the other hand, the topology of the 
network, its interconnection and damping laws at the vertices, and the rational 
dependence of the transit times on the network edges [1,5,8,10,29,30]. A case 
of special interest is when some coefficients of the system are time-dependent and 
switch arbitrarily within a given set [3, 15,24]. 

In this paper, we address stability issues first for transport systems with time- 
dependent transmission conditions and then for wave propagation on networks with 
time-dependent damping terms. When the time-dependent coefficients switch ar¬ 
bitrarily in a given bounded set, we prove that the stability is robust with respect 
to variations of the lengths of the edges of the network preserving their rational 
dependence structure (see Corollary 4.11 for transport and Corollary 5.15 for wave 
propagation). Such robustness enables us to get the main result of the paper, namely 
a necessary and sufficient criterion for exponential stability of wave propagation on 
networks: exponential stability holds for a network if and only if it is a tree and 
the damping is bounded away from zero at all external vertices but at most one 
(Theorem 5.16). 

We address these issues by formulating them within the framework of non- 
autonomous linear difference equations 

N 

u{t) = '^Aj{t)u{t- Aj), «(!) e C'^, (Ai,..., Ajv) G (K+)^. (1) 

1=1 

This standard approach relies on the d’Alembert decomposition and classical trans¬ 
formations of hyperbolic systems of PDEs into delay differential-difference equa¬ 
tions [6,9,12,19,22,27]. Here, stability is meant uniformly with respect to the 
matrix-valued function A(-) = (Ai(-),... ,An{-)) belonging to a given class A. 

In the autonomous case. Equation (1) has a long history and its stability is com¬ 
pletely characterized using Laplace transform techniques by the celebrated Hale- 
Silkowski criterion (see e.g. [4, Theorem 5.2], [16, Chapter 9, Theorem 6.1]). The 
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latter can be formulated as follows: if Ai,..., Ajv are rationally independent, then 
all solutions of u{t) = ^ju{t — Kj) tend exponentially to zero as t tends to 

infinity if and only if /Ohs(^) < 1) where 

(2) 

and p(-) denotes the spectral radius. Notice that the latter condition does not 
depend on Ai,..., Ajy. The Hale-Silkowski criterion actually says more, namely the 
striking fact that exponential stability for some Ai,..., Atv rationally independent 
is equivalent to exponential stability for any choice of positive delays Li,..., Ljy. 
This criterion can be used to evaluate the maximal Lyapunov exponent associated 
with u{t) = — Aj), i.e., the infimum over the exponential bounds for 

the corresponding semigroup. A remarkable feature of the Hale-Silkowski criterion 
is that, contrarily to the maximal Lyapunov exponent, it does not involve taking 
limits as time tends to infinity. An extension of these results has been obtained 
in [21] for the case where Ai,..., Ajv are not assumed to be rationally independent. 

The non-autonomous case turns out to be more difficult since one does not have 
a general characterization of the exponential stability of (1) not involving limits 
as time tends to infinity. To illustrate that, consider the simple case N = 1 of 
a single delay and A = f8) where *8 is a bounded set oi d x d matrices. 

Then the stability of (1) is equivalent to that of the discrete-time switched system 
Un+i = AnUn where An € and it is characterized by the joint spectral radius of 

(see for instance [18, Section 2.2] and references therein) for which there is not 
yet a general characterization not involving limits as n tends to infinity. 

Up to our knowledge, the only results regarding the stability of non-autonomous 
difference equations were obtained in [23], where sufficient conditions for stability 
are deduced from Perron-Frobenius Theorem. Our approach is rather based on a 
trajectory analysis relying on a suitable representation for solutions of (1), which 
expresses the solution u{t) at time t as a linear combination of the initial condition 
uo evaluated at finitely many points identified explicitly. The matrix coefficients, 
denoted by 0, are obtained in terms of the functions Ai (•),..., AAr(-) and take 
into account the rational dependence structure of Ai,..., Ajv (Proposition 3.14). 
This representation provides a correspondence between the asymptotic behavior of 
solutions of (1), uniformly with respect to the initial condition uq and A(-) G A, 
and that of the matrix coefficients 0 uniformly with respect to A(') £ A (Theorem 
3.21). In the case where A = L°“(R, *8) for some bounded set *8 of A^-tuples of dx d 
matrices, we extend the results of [21], replacing pus in the Hale-Silkowski criterion 
by a generalization p, of the joint spectral radius. As a consequence of our analysis, 
and despite the lack of a closed and delay-independent formula for p analogous 
to (2), we are able to show that stability for some N-tuple A = (Ai,...,Ajv) 
is equivalent to stability for any choice of A^-tuple (Li,... jLat) having the same 
rational dependence structure as A (Corollaries 3.29 and 3.35). 

The structure of the paper goes as follows. Section 2 provides the main nota¬ 
tions and definitions used in this paper. Difference equations of the form (1) are 
discussed in Section 3. We start by establishing the well-posedness of the Cauchy 
problem and a representation formula for solutions in Sections 3.1 and 3.2. Stability 
criteria are given in Section 3.3 in terms of convergence of the coefficients and spec¬ 
ified to the cases of shift-invariant classes A and arbitrary switching. In the latter 
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PHS(^) = 


max 


(@i,...,0A^)G[O,27r]^ 




,6 
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case, we provide the above discussed generalization of the Hale-Silkowski criterion. 
Applications to transport equations are developed in Section 4 by exhibiting a cor¬ 
respondence with difference equations of the type (1). Thanks to the d’Alembert 
decomposition, results for transport equations are transposed to wave propagation 
on networks in Section 5. The topological characterization of exponential stability 
is given in Section 5.3. 

2. Notations and definitions. In this paper, we denote by N and N* the set of 
nonnegative and positive integers respectively, K+ = [0,-|-oo) and = (0,-|-oo). 
For a, 6 G K, let |o, 6] = [a, b] fl Z, with the convention that [a, 6] = 0 if a > 6. The 
closure of a set F is denoted by T. If a: G R and F C R, we use a; -I- T to denote 
the set {x + y \ y G F}. 

We use and 6ij to denote, respectively, the cardinality of a set F and the 
Kronecker symbol of i,j. For a: G K, we use x± to denote max(±a::, 0) and we 
extend this notation componentwise to vectors. For x G we use a^min and x„iax 
to denote the smallest and the largest components of x, respectively. 

If AT is a subset of C and d,m G N, the set of d x m matrices with coefficients in 
K is denoted by A4d,miK), or simply by A4d{K) when d = m. The identity matrix 
in Md{C) is denoted by Id^. We use ei,..., to denote the canonical basis of 
i.e., Ci = {5ij)jeii^d} for i G |l,d]]. For p G [l,-|-oo], l-l^ indicates both the f^-norm 
in and the induced matrix norm in A4d{C). We use p{A) to denote the spectral 
radius of a matrix A G A4d(C), i.e., the maximum of |A| with A eigenvalue of A. 
The range and kernel of a matrix A are denoted by Ran A and Ker A respectively, 
and rk(A) denotes the dimension of Ran A. Given Ai,..., A^r G Md{C), we denote 
by nti ^ i the ordered product Ai ■ ■ ■ Aj^. 

All Banach and Hilbert spaces are supposed to be complex. For p G [1, -foo], we 
use to denote the usual Lebesgue spaces of p-integrable functions and the 
Sobolev spaces of A:-times weakly differentiable functions with derivatives in . 

A subset A of L“j.(R, Ald(C)^) is said to be uniformly locally bounded if, for 
every compact time interval I, sup^g _4 || Ajj^oo^j. Md{C)^) finite. We say that A is 
shift-invariant if A(- 1) G A for every A G A and < G R. 

Throughout the paper, we will use the indices S, r and w in the notations of 
systems and functional spaces when dealing, respectively, with difference equations, 
transport systems and wave propagation. 

3. Difference equations. Let N, d G N*, A = {Ai,..., An) G (R^)^, and con¬ 
sider the system of time-dependent difference equations 

N 

Es{A,A)-. u(t) =^Aj(t)u(f-Aj). (3) 

4 = 1 

Here, u(t) G and A = (Ai,..., A^v) : R —>■ A4d(C)^. 

3.1. Well-posedness of the Cauchy problem. In this section, we show existence 
and uniqueness of solutions of the Cauchy problem associated with (3). We also 
consider the regularity of these solutions in terms of the initial condition and A(-). 

Definition 3.1. Let Uq : [— Amax, 0) —)> and A = (Ai,..., An) ■ R —>■ Md{'C)^ . 
We say that u : [—Amaxi +oo) —>■ C'’* is a solution of S, 5 (A, A) with initial condition 
Uq if it satisfies (3) for every t G R+ and u{t) = uo{t) for t G [—A^ax, 0). In this 
case, we set, for t>0,ut = u{- + t)|[-A„ax.o)- 
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We have the following result. 

Proposition 3.2. Let uq : [—AmaxjO) —>■ and A = (Ai,..., Ajv) : K. —>■ 

Aid{C)^. Then Tis{K^A) admits a unique solution u : [—Amax)+oo) —>■ with 

initial condition uq. 

Proof. It suffices to build the solution u on [—Amax, Amin) and then complete its 
construction on [Aminj +oo) by a standard inductive argument. 

Suppose that u : [—Amax,Amin) —is a solution of ^^(A, A) with initial 
condition uq. Then, by (3), we have 

f "" 

I ^ ^ Aj{t)uo{t — Aj), if 0 < t < Amin, 

Up) = ] 

[ Uo{t), if - Amax < t < 0. 

Since the right-hand side is uniquely defined in terms of Ug and A, we obtain the 
uniqueness of the solution. Conversely, if u : [—Amax, Amin) —>■ is defined by (4), 
then (3) clearly holds for t G [—Amax, Amin) and thus u is a solution of ^^(A, A). □ 

Definition 3.3. For p € [l,-|-oo], we use to denote the Banach space = 
LP([—Amax, 0], C^^) endowed with the usual L^-norm denoted by [[-[[p. 

Remark 3.4. If ug, vg : [—Amax, 0) —>■ are such that ug = vg almost everywhere 
on [—Amax, 0) and A,B :!&—>■ A4d(C)^ are such that A = B almost everywhere on 
K+, then it follows from (4) that the solutions u, v : [—Amax, +oo) —>■ C‘^ associated 
respectively with A, ug and B, vg satisfy u = v almost everywhere on [—Amax, +oo). 
In particular, for initial conditions in X^, p G [l,-|-oo], we still have existence and 
uniqueness of solutions, now in the sense of functions defined almost everywhere. If 
moreover A G L^^(R, A4d(C)^), it easily follows from (4) that the corresponding 
solution u of ^^(A, A) satisfies u G Af^cd^^max, +oo), C^). 

Proposition 3.5. Let p G [l,-|-oo), ug G Xp, A G L^^{R, MdiC)^), and u he the 
solution of'Es{A, A) with initial condition ug. Then the X^-valued mapping 1ut 
defined on ]R_|_ is continuous. 

Proof. By Remark 3.4, ut G Xp for every t G K+. Since ut{s) = u{s + t) for 
s G [—Amax, 0), the continuity of 1 1 —>■ rt* follows from the continuity of translations 
in LP (see, for instance, [25, Theorem 9.5]). □ 

Remark 3.6. The conclusion of Proposition 3.5 does not hold for p = -foo since 
translations in L°° are not continuous. 

3.2. Representation formula for the solution. When t G [0,Amin), Equation 
(4) yields u(t) in terms of the initial condition ug. If t > Amin, we use (3) to express 
the solution u at time t in terms of its values on previous times t — Aj, and, for each 
j such that t > Aj, we can reapply (3) at the time t — Aj to obtain the expression 
of u{t — Aj) in terms of u evaluated at previous times. By proceeding inductively, 
we can obtain an explicit expression for u in terms of ug. For that purpose, let us 
introduce some notations. 

Definition 3.7. i. An increasing path (in N^) is a finite sequence of points 

(Qfc)fe=i such that, for k G |1, n—1|, q/c+i is obtained from q/j by adding 
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1 to exactly one of the coordinates of q^. For n G N* and v = {vi,... ,Vn) G 
|1, iV]"', we use to denote the increasing path defined by 


fe-i 

Pvik) = 

i=i 

ii. For n G \ {0}, we use Vn to denote the set 


V„ = 


e p«(|n|i + l) =n|, 


i.e., Vn can be seen as the set of all increasing paths from 0 to n. 
iii. For A = (Ai,...,Ajv) :R^ Md(C)^, A = (Ai,...,Ajv) G (K;)^, n g 
and t G M, we define the matrix G Afd(C) inductively by 

fO, ifnGZ^\N^, 

Idd, if n = 0, 

ifnGN^\{0}. 


-^A,A 


= 


(5) 


fc=i 


We omit A, A or both from the notation when they are clear from the 
context. 


The following result provides a way to write as a sum over Vn and as an 
alternative recursion formula. 


Proposition 3.8. For every n G N'^ \ {0} and t gR, we have 


and 


-n.t = ri ^ • P«(^)) 

V^Vn k — 1 


^A,A 

^n,t 


N 

- A • n + Afe). 




( 6 ) 

(7) 


Proof. We prove (6) by induction over |n|^. If n = for some i G |1, A^], we have 

1 

k — 1 


Let i? G N* be such that (6) holds for every n G with |n|j^ = R and t G K. If 
n G is such that |n|j^ = i? + 1 and t G M, we have, by (5) and the induction 
hypothesis, that 


^n,t 


N 

Ak{t)^n-ek,t-Ak 


fc=l 

nfc>l 


E E 


1— 1 Pe Vn — e 
>1 


rik 


l“li 

v^Vn r=l 


|n|i-l 

Ak{t) p (t - Afc - A • p„(r)) 


where we use that Vh = {{k,v) \ k G |1,A^], nt > 1, v G Vn-ei,} and that 
Gfc + P«(r) = P(k,v){^ + !)• This establishes (6). 
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We now turn to the proof of (7). Since ^ej,t = (7) is satisfied for n = e^-, 

j G For n G with |n|^ > 2, the set Vn can be written as 

Vn = {{v,k) I k G |l,iV]|, rife > 1, n G Vh-et}, 
and thus, by (6), we have 

N 


|n|l-l 

Ay^ (t - A ■ Py{r)) 

r—1 


m=E e 

k—1 V^Vn- 
nfc>l 

N 

= E E 

k—1 V^Vn- 
nfc>l 

N 

— ^ ^ ^n—A ■ n + A/j;)- 


(f-A-p„(r)) 

r—1 


Akit- A- p„(|n|J) 


Ak{t - A ■ n + Ak) 


□ 


fc=i 


In order to take into account the relations of rational dependence of Ai,..., Ajv G 
in the representation formula for the solution of ^^(A, A), we set 

Z(A) = {nGZ^ |A-n = 0}, 

F(A) = {L G I Z{A) c Z{L)}, V+{A) = V{A) n (R^)^, 

W{A) = {L G R^ I Z{A) = Z{L)}, W+{A) = W{A) n (R;;)^. 

Notice that IF(A) c y(A) and IF(A) = {L G V{A) \ V{L) = F(A)}. 

The point of view of this paper is to prescribe A = (Ai,..., Ajv) G (R!j_)^ and 
to describe the rational dependence structure of its components through the sets 
Z{A), F(A), and W{A). Another possible viewpoint, which is the one used for 
instance in [21], is to fix B G Ad N,h{^) consider the delays A = (Ai,..., An) G 

Ran_B fl (R^)^. We show in the next proposition that the two points of view are 
equivalent. 

Proposition 3.9. Let A = (Ai,...,A 7 v) G (R!]_)'^. There exist h G |1,A^], £ = 
{£i,... ,£h) G (R^)^ with rationally independent components, and B G A4N,h(^) 
with rk(il) = h such that A = B£. Moreover, for every B as before, one has 

V(A) = RanR, 

/o') 

W (A) = {B{£'^,... ,£'f,) \ £[,..., £'f,^ are rationally independent}. 

In particular, IF(A) is dense and of full measure in F(A). 


Proof. Let V = Span^lAi,..., A^v}, h = dim^V, and {Ai,...,Aft,} be a basis of 
V with positive elements, so that A = Au for some A = {oij) G M-N.hiQ) and 
u = (Ai,...,A;j) G (R^)^. For G R^ \ {0}, we denote by Py the orthogonal 
projection in the direction of v, i.e., Py = vv"^/{vl^. 

Since is dense in R^, there exists a sequence of vectors ..., rh^n) 

in (Q!j_)^ converging to n as n —>■ +oo, and we can further assume that the sequence 
is chosen in such a way that |P„„ — Pu \2 ^ 1/n^ for every n G N*. 

For n G N*, we define Tn = i {Idh—Pun)- This operator is invertible, 

with inverse T~^ = Py^ + n (Id/j —Py^). Furthermore, both Ty and T~^ belong to 
MhiQ). For i G |l,h|, we have 

(T~^ei)'^u = eJPy^u + nej{ldh -Py„)u = eJPy^u + nej{Py - 
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and thus (T„ —>■ eju = Ai as n —>■ +oo. Since Ai, ... ,Xh > 0, there exists 

no G N* such that 

{T~^ei)'^u >0, Vi G |1, /i|, Vn > uq. (9) 

For i G |1,-/V], let ai = {aij)jeii^hj G For each i G |1,-/V]|, we construct the 
sequence in by setting ai^n = Tntti. It follows from the definition of 

T„ that ai^n converges to PuOi = “1^ as n —>■ +oo. Since u^ai = = 

Ai > 0 and the components of u are positive, we conclude that there exists ni > no 
such that ai,ni G (Q+)^ for every i G |1, Af|- 

Let £ = {T~^)'^u. By (9), £i = {T~^ei)'^u > 0 for every i G Since the 

components of u are rationally independent, £i,... ,£h are also rationally indepen¬ 
dent. Let bij G Q+, i G |1,A^|, j G be such that ai^m = • Hence, 

for i G [1, A^l, 

h h 

U Oii 'll Q^z,ni ^ ^ ^ij'^ '^ni ^ ^ ^^3^3' 

We then get the required result up to multiplying B = (6y ) by a large integer and 
modifying £ in accordance. 

We next prove that (8) holds for every B as before. (Our proof actually holds 
for every B G AiN,h{Q) with rk(i?) = h such that A = B£ for some I G (K.!j_)^ 
with rationally independent components.) First notice that Z{K) = {n G I n € 
KerS"^}. Indeed, n G Z{K) if and only if n is perpendicular in to B£^ which 
is equivalent to B = 0 since £i,... ,ih are rationally independent. Moreover, re¬ 
mark that Keri?"’’ = (RanB)-'* admits a basis with integer coefficients since RanB 
admits such a basis. To see that, it is enough to complete any basis of Rani? in 
hy N — h vectors in and find a basis of (Rani?)-^ by Gram-Schmidt orthogo- 
nalization. We finally deduce that SpanR(Z(A)) = (Rani?)-^. Since by definition 
R(A) = Z(A)-*-, we conclude that V{A) = Rani?. As regards the characteriza¬ 
tion of W{A), an argument goes as follows. Let L G R(A), so that L = B£' for 
a certain £' G The components of £' are rationally dependent if and only if 
dimQSpanQ{Li,...,L tv} < h, i.e., dimRR(L) < dimRR(A), which holds if and 
only if L ^ W{A). □ 

We introduce the following additional definitions. 

Definition 3.10. Let A = (Ai,...,A 7 v) G We partition |1,A^] and 

according to the equivalence relations ^ and ~ defined as follows: * ^ j if A^ = Aj 


and n Ri n' if A • n = A • n'. We use [•] to denote equivalence classes of both ^ 
r; and we set J = |1, A^]/ ~ and Z = j r;. 

For A : R — MdiC)^ , L G V+(A), [n] G Z, [i] G J, and t G R, we define 

and 

^L,A,A _ 
“[n],t ~ 

n'G[n] iG[i] 

(10) 

and 

pvL,A,A 


(11) 


bleJ 


Remark 3.11. The expression for given in (10) is well-defined since, thanks 

to (5), all terms in the sum are equal to zero except finitely many. The expression 
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for given in (11) is also well-defined since, for every L G V+(A), it i ^ j and 

n Ri n', one has L,; = L, and L ■ n = L ■ n'. In addition, notice that ^ 0 


only if [n] fl is nonempty, and, similarly, 7 ^ 0 only if [n] fl (N^ \ { 0 }) 

is nonempty. Another consequence of the above fact and (11) is that 7 ^ 0 

only if f > 0, since [n — ej] fl = 0 whenever [n] € Z and [j] G J are such that 
L ■ n — Lj <0. 

Notice, moreover, that the matrices S, A and 0 depend on A only through Z{A). 
Hence, if A' G IH+(A) (i.e., Z(A') = Z(A)), then 

lA _ AA' / 


qn],t 


SL,A,A _ 5L,A',A 
“ “fnl.t > 




qL A A ^ 0T,A',A_ 
[n],t [n],t 


From now on, we fix A = (Ai,..., Ajv) G (M^)^ 


and our goal consists of deriving 
a suitable representation for the solutions of 'Zs{L,A) for every L G F+(A). Even 
though the above definitions depend on A, L G V+(A) and A, we will sometimes 
omit (part of) this dependence from the notations when there is no risk of confusion. 
Let us now provide further expressions for ^ . 

Proposition 3.12. For every L G 1 + (A), A : R —7 Md{C)^, n G \ {0}, and 
< G R., we have 


“[n],t ~ ^ b] ’ 

m^j 

and 


qn],t 


= E 

lj]ej 


‘[n-ej],t b] 


At] {t — L ■ n + Lj), 


E E (12) 

[n]nN-^ pG ^—1 

Proof. We have, by Definition 3.10 and Equation (5), that 

N N 

“Mg = E ““'G “ E E“ E E 


n'G[n] 

N 


n'G[n] j=l 


1=1 


n'G[n] 


N 




1=1 


1=1 


— E E A(t)‘=‘[n-ei].t-Li — E I E A(t) I ^[n-ej],t-Lj 

[l]etriG[i] [1]G.7 \ie[i] 

= E ■ 

blej 

The second expression is obtained similarly from Definition 3.10 and Equation (7) 
and the last one follows immediately from (6) and (10). □ 

Let us give a first representation for solutions of 'Zs{L, A). 


Lemma 3.13. Let L G 


p* \N 


A = (Ai,...,A^) : 


AiddC)^, and uq : 


[—Lniax, 0) —I C“. The eorresponding solution u : [—Lmax,+oo) —t C“ of Tis{L, A) 


is given for t > 0 by 

uit) = 


E 


-n-e^,tAj(t -L-n + Lj)uo(t-L-n). 


(nj)eN'^x[l.Afl 
— Z/j<t—L-n<0 


(13) 
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Proof. Let u : [—Lmax, +c») —>■ be given for t > 0 by (13) and u{f) = uo{t) for 

t £ [—LmaxjO). Fix t>0 and notice that 

N 

Lj) 

1=1 

N 

= Y1 H Lj- L-n + Lk)uo{t-Lj- L-n) 

(n,fc)GN^x|l,Ar| 

—Lfc<t—Lj —L-n<0 

nfc>l 

N 

+ '^ Lj). (14) 

1=1 

t<Lj 

Consider the sets 

Bi{t) = {(n, fc, j) e X |l,iVp 1 1 > Lj, -Lk <t - Lj - L ■ n < 0, nfc > 1}, 
B2{t) = {jG 11, Nj\t<L,}, 


Ci{t) = {{n,k,j) G N^x|l, \-Lk < t-L n < 0, > 1, Uj > l+Sjk, n ^ Cfe}, 


C'2(t) = {(n, fc) G X |1,1V| I - Lfe < t - L • n < 0, n = Ck }, 
and the functions ipi : Bi{f) —1 Ci{t), i G {1,2}, given by 


ipi{n,k,j) = {n + ej,k,j), ip^ij) = (e^.j). 


One can check that pi and (p 2 are well-defined and injective. We claim that they 
are also bijective. For the surjectivity of pi, we take {n,k,j) G C'i(t) and set 
m = n — Cj . Since nj > 1, one has m G N'^. Since n/c > 1, rij > 1 -f Sjk, one has 
t > L ■ n — Lk > Lj + Lk — Lk = Lj. The inequalities —Lk <t — Lj — Lm<Q 
and rife > 1 are trivially satisfied, and thus (m, fc,j) G Bi{t), which shows the 
surjectivity of ipi since one clearly has (/^i(m, k,j) = (n, k,j). For the surjectivity 
of (/? 2 , we take (n, k) G C 2 (t), which then satisfies n = and t < L ■ n = Lk. This 
shows that k G B 2 {t) and, since p> 2 {k) = (n, k), we obtain that ip 2 is surjective. 

Thanks to the bijections pi, (p 2 , and (5), (14) becomes 

N 

'^Aj{t)u{t- Lj) 

1=1 

N 

= H Aj{t)E^-:^^^_^. ^_^Ak{t-L-n + Lk)uo{t-L-n) 

(n,/c)GN'^xfl,Afl 1=1 

-Lfc<t-L.n<0 nj>l+Sjk 
nfc>l, n/efc 

+ ^ Akit - L ■ n + Lk)uo{t - L ■ n) 

(n,fc)GN'^x[l,Af] 

— Lk<t—L-n<0, 
n—ek 


STABILITY OF DIFFERENCE EQUATIONS AND APPLICATIONS 


11 


= ^n’^^^^tAk{t-L-n + Lk)uo{t-L-n) 

{n,k)eN^ xll,N} 

— Lfc<t —L-n<0 
»fe>l, n/efc 

+ ^ E^^f^Ak{t-L-n +Lk)uo{t-L-n) 

(n,fe)GN'^x[l,Afl 
— Lfc<t —L-n<0, 

n=efc 

= ^t-ek,t^k{t-L-n + Lk)uo{t-L-n) = u{t), 

(n.fc)GN"x[l,Afl 
— Lfc<i—L-n<0 

which shows that u satisfies (3). □ 

We can now give the main result of this section. 

Proposition 3.14. Let A e L G V+{A), a : K ^ Md(C)^, and uo : 

[—Amax, 0) —>■ C'^. The corresponding solution u : [—Lmax,+oo) —t of Yjs{L^ A) 
is given for t > 0 by 

u{t)= L-n), (15) 

[n]ez 

t <C Z/*ri*Ct-|-Lmfl,x 

where the coefficients 0 are defined in (11). 

Proof. Equation (15) follows immediately from (13) and from the fact that the 
function ip : x |1, Af| —>• Z x x J x |1, A^] given by ip{n,j) = ([n], n, [j], j) is 

a bijective map from {(n, j) G x |1, N} \ —Lj < t — L-n< 0} to {([m], n, [i],j) G 
Z X X J" X |1, A^JI I n G [m], j G [i], t < L • n < t + L„iax, L ■ n — Lj < t} for 
every t G M. □ 

3.3. Asymptotic behavior of solutions in terms of the coefficients. Let us 

fix a matrix norm I-] on Ald(C) in the whole section. Let ( 71,(72 > 0 be such that 

(7i |7l|p < |A| < C 2 |A|p, VAg Ald(C), VpG [l,+oo]. (16) 

Let A be a uniformly locally bounded subset of L“j,(R, A4d{C)^). The family of 
all systems Tis{L, A) for A G A is denoted by (L, A). We wish to characterize the 
asymptotic behavior of ^^(L, A) (i.e., uniformly with respect to A G A) in terms 
of the behavior of the coefficients ^ and 0 [n],t- Foi' that purpose, we introduce 
the following definitions. 

Definition 3.15. Let L G (R!^)^. 

i. For p G [l,+oo], we say that 'Es{L,A) is of exponential type 7 G K in if, 
for every e > 0, there exists A > 0 such that, for every A G A and uq G X^, 
the corresponding solution u of 'Es{L,A) satisfies, for every t > 0 , 

We say that ^^(L, A) is exponentially stable in Xp if it is of negative exponen¬ 
tial type. 

ii. Let A G (R!J_)^ be such that L G V+(A). We say that T,s{L,A) is of (0, A)- 
exponential type 7 G K if, for every e > 0, there exists K > 0 such that, for 
every A G A, n G N-^, and almost every t G (L • n — Lmax, A ■ n), we have 


(^L,K,A 
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iii. Let A s (R^)^ be such that L G V+(A). We say that T,s[L,A) is of (S, A)- 
exponential type 7 G R if, for every e > 0, there exists K > 0 such that, for 
every A G .4, n G and almost every t G R, we have 


"L,K,A 

“[n],t 




iv. For p G [1, +c»], the maximal Lyapunov exponent of T,s{L, A) in is defined 


as 


Xp{L,A) = limsup sup sup 


InllwtlL 


*->■+00 AeA UQgXp 
lkolL=i 


where u denotes the solution of Yis{L, A) with initial condition uq. 


Remark 3.16. Let L G (Rl^)^ and /r G R. For every A : R — >■ Ald(C)^ and u 
solution of Tjs{L, A), it follows from (3) that 1 1 —>■ e^*u{t) is a solution of the system 
X^siL, {e^^^Ai ,..., Aim)). As a consequence, if A C L“^(R,44^(0)'^) and 

Ap. = {(e'^-^^Ai,..., e^^^Aw) | A = (Ai,..., A^) G A}, 

one has Ap(L, A^) = Ap(L, A) + p. 


The link between exponential type and maximal Lyapunov exponent of Tis{L, A) 
is provided by the following proposition. 

Proposition 3.17. Let L G (RlJ.)'^, A be uniformly locally bounded, and p G 
[l,+oo]. Then 

Xp{L,A) = inf {7 G R | A) is of exponential type 7 in Xp}. 

In particular, A) is exponentially stable if and only if Xp{L,A) < 0. 

Proof. Let 7 G R be such that ^^(L, A) is of exponential type 7 in X^. It is clear 
from the definition that Xp{L,A) < 7 . We are left to prove that Iis{L,A) is of 
exponential type Xp{L,A) when the latter is finite. Let e > 0. From the definition 
of Ap(L,A), there exists to > 0 such that, for every t > to, A G A, and uq G Xp, 
one has 

Since A is uniformly locally bounded, by using (13) and ( 6 ), one deduces that 
there exists K > 0 such that, for every t G [0,to], A G A, and ug G Xp, one has 
IKtIlp < at II Mo lip. Hence the conclusion. □ 


Remark 3.18. Similarly, one proves that, for A G (R1J_)^ and L G F+(A), 


lim sup sup ess sup 

L/'ti —^+00 A^.A. 


In 

qL,A,A 


[n],i 


t 


= inf {7 G R I T,s{L,A) is of (0, A)-exponential type 7 } 


and 


lim sup sup ess sup 

L'ti —^-|-cxD A^A. 


In 

"L,A,A 


[n],t 


L • n 


= inf {7 G R I Iis{L,A) is of (^, A)-exponential type 7 }. 
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3.3.1. General case. The following result, which is a generalization of [ 8 , Proposition 
4.1], uses the representation formula (15) for the solutions of Y,s{L,A) in order to 
provide upper bounds on their growth. 

Proposition 3.19. Let L G P+(A). Suppose that there exists a continuous func¬ 
tion /:]&—>■ such that, for every A G A, n G , and almost every 
t G (L • n — Tmax) L ■ n), one has 



(17) 


Then there exists C > 0 such that, for every A G A, p G [1, +oo], and uq GX^, the 
corresponding solution u of A) satisfies, for every t>Q, 


llwtllp < (^(t + 1 )^ ^ max /(s)||Mo||p- 

SG [t i^max 


(18) 


Proof. Let A G A, p G [l,+oo), uq G X^, and u be the solution of T,s{L,A) 
with initial condition uq. For t G M+, we write Jtf = {M G .Z | t < L • n < 
t + Lmax, [n] n ^ 0} and Vt = #yt. Thanks to Proposition 3.14, Remark 3.11, 
and (17), we have, for t > Tmax, 



[n]Gj^s 



We clearly have Yt < #{n G \t<L-n<t-\- Tmax}- For n G N^, we denote 
Cln = {a; G \ Hi < Xi < Hi I for every i G |1, A^]}. This defines a family of 
pairwise disjoint open hypercubes of unit volume. Thus 


f \ 



< Vol{x G (R+)^ I f < L • a; < t + |L|^ + Lmax}- 


N-l 


Then there exists C 3 > 0 only depending on L and N such that Yt < C' 3 (t + 1) 
Thus, 
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max f{sY 

•SG [i — -£/max,i] 



|Mo(s)|pds. 


Similarly, there exists (74 > 0 only depending on L and N such that, for every 
t G M+ and s G [-Lmax, 0], Yt-L^^^-s < C^t + 1)^“\ yielding (18) for t > Lmax- 
One can easily show that, for 0 < t < Lmax, we have ||ut||p < C ||ito||p for some 
constant C' independent of p and uq, and so (18) holds for every t > 0. The case 
p = +00 is treated by similar arguments. □ 


When L G W7|-(A), we also have the following lower bound for solutions of 
YsiL,A). 


Proposition 3.20. Let L G W+(A) and /:!&—>■ be a continuous function. 
Suppose that there exist A G 7l, no G N^, and a set of positive measure S C 
(L ■ no — Amax, A • no) such that, for every s € S, 


0 


L,A,A 

[no],s 


> /(«)■ 


(19) 


Then there exist a constant (7 > 0 independent of f, an initial condition ug G 
A°“([— Aniax, 0],C*^), and t > 0, such that, for every p G [l,+oo], the solution u of 
Tis{L,A) with initial condition uq satisfies 


\Ut 


\ >C min 

sC [t Lma. 




Proof. According to Remark 3.11, one has for every [n] G Z and 

s G M, and therefore we assume for the rest of the argument that A = L and we 
drop the upper index A, A, A. 

For s G S', one has |0[no],s|oo > Clf^fis), where C 2 is defined in (16). Using (19) 
and Remark 3.11, one derives that S C [0, + 00 ). 

For every s G S, one has 

d 

C 2 /(s) < |0[no],s|^ < ^ , 

i=i 

and thus there exist jo G |l,(il and a subset S C S of positive measure such that, 
for every s G S and p G [1, + 00 ], one has 

CY^d-^f{s) < |0[no].sejo|^ < |0[no].aejo|p. (20) 

In order to simplify the notations in the sequel, we write S instead of S. 

Let to G S \ {0} be such that, for every e > 0, (to — £,to + e) fl S has positive 
measure. Let <5 > 0 be such that 


f 1 

2i5< min< 2to,A-no-to,to-A-no + Amax, mm |A-(n-no)|>. 

nGN" 

L.(n-no)/0 

Such a choice is possible since to G (A • no — Amax, A • no), tp G S' \ {0} C K+, and 
{A • n I n G N^} is locally finite. 

Let Si = (S — to) n {—S,S), which is, by construction, of positive measure, and 
/i : R —>■ R be any non-zero bounded measurable function which is zero outside Si. 
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Define uq : [-Lmax,0) by 

uq(s) = n{s-to+ L ■ no)ejo 

and let u be the solution of Tis{L,A) with initial condition uq. For s G (—(5,(5), we 
have to + s > 0 since to > <5- By Proposition 3.14, one has 

u(to + s)= ^ 0[„],to+s/r(s + L-(no-n))ejo. (21) 

[n]G2 

tQ-\-S<iL'll'^tQ-\-S-\-Lxxia.x 

If L • n ^ L • no, we have |L • (n — no)| > 2d, and so |s + L • (no — n)| > S, which 
shows that /r(s + L ■ (no — n)) = 0. Hence, Equation (21) reduces to u(to + s) = 
0 [no],to+st^(®) 6 io- We finally obtain, using (20) and letting t = to + d, that, for 
p G [l,+oo), 


Ltllp > ||MtollL([-M],C‘^) ^ / \u{to + s)\lds= f |0[no].to+seio|p L(s)|^ds 

^ Si J Si 

>C^Pd-P [ f{to + sr\Ksrds>C^Pd-P min /(s)^ ||uo||^ . 

Js, se[t-imax,t] 


( 22 ) 


A similar estimate holds in the case p = +oo, which concludes the proof of the 
proposition. □ 


As a corollary of Propositions 3.19 and 3.20, by taking / of the type /(t) = 
one obtains the following theorem. The last equality follows from Propo¬ 
sition 3.17 and Remark 3.18. 


Theorem 3.21. Let A G (K.^)^ and A be uniformly locally bounded. For every 
L G 1+(A), if 'Fs{L, A) is of (0, K)-exponential type 7 then, for every p G [1, -boo], 
it is of exponential type 7 in X^. Conversely, for every L G W+(A), if there exists 
p G [l,+oo] such that Iis{L,A) is of exponential type 7 in X^, then it is of (0, A)- 
exponential type 7 . Finally, for every L G W+(A) and p G [1, -boo], 


Xp{L,A) = limsup sup ess sup 

L-n —>^-|-oo AG-A —Z/maxjL-n) 


In 

0L A A 


[n],t 


t 


(23) 


Remark 3.22. It also follows from Proposition 3.19 that, in the first part of the 
theorem, the constant AT > 0 in the definition of exponential type of 'FsiL,A) can 
be chosen independently of p G [1, -boo]. Moreover, the left-hand side of (23) does 
not depend on p and its right-hand side does not depend on A. 


3.3.2. Shift-invariant classes. We start this section by the following technical result. 

Lemma 3.23. For every A G (K^)^, L G V+(A), A : R —^ , n G and 

t, r G K, we have 

^L,A _^L,A(.+r) . SL.A.A _ eL.A,A(.+r) 

^n.t+T — ^n,t “[n],t-|-T “ “[n],t 

Proof. The first part holds trivially if n G \ or if n = 0, for, in these cases, 
it follows from (5) that does not depend on t and A. If n G \ {0}, the 
conclusion follows as a consequence of the explicit formula ( 6 ) for The second 

part is a consequence of the first and ( 10 ). □ 
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We next provide a proposition establishing a relation between the behavior of 
S[n],t and 0[n],t- Notice that, if a subset A of L^^(R, A4d(C)^) is shift-invariant, 
then A is uniformly locally bounded if and only if it is bounded. 


Proposition 3.24. Let A he a bounded shift-invariant subset of L°°(M.,A4d{C)^), 
L G V+(A), and /:!&—>■ be a continuous function. Then the following assertions 
hold. 


i- If 


0 


L,A,A ^ holds for every A G A, n G , and almost every t G 
(L • n — Lmax, L ■ n), then, for every A G A, n G \ {0}, and almost every 


t G M, one has 


'[n],t 


ii. If 


^L,A,A 

“[n],t 


< /(L • n) holds for every A G A, n G , and almost ev¬ 


ery t G R, then there exists a constant C > 0 such that, for every A G A, 


n G N-^, and almost every t G (L ■ n — Lmax, L ■ n), one has 
f{s). 


0 


L,A,A 

n] 


< 


Proof. We start by showing (i). Let A G A and n G \ {0}. For every k G "L, 
there exists a set iVfc C [L • n — Lmin, L • n) of measure zero such that, for every 
t G [L • Yl L^nin: A • n) \ 


fcAmin) 

= 

E -L-n + L,) 


= 

qZ/,A,A(- A:Z;min) 

< fit), 


where we use Proposition 3.12, the fact that L ■ n — Lj < L ■ n — Lmin < t for every 
[i] G J^ and Equation (11). 

Let N = UfcGz(-^fe ~ kLmin), which is of measure zero. For t G R \ N', let 
fc G Z be such that t G [L • n — (fc -|- l)Lmin, L ■ n — fcLmin), so that t -\- kL^in G 
[L • n — Lmin, L ■ n). Since t ^ N, we have t -|- fcLmin ^ Nk, and so, using Lemma 
3.23, we obtain that 


^L,A,A 


i^Z/,A,A(- AlZ^niin) 

“[n],t 


^[n],£ + fcLmin 


< f{t-\-kLn,in) < max f{s). 

sG[L-n-Lmln.Z-n] 


Let us now show (ii) . Without loss of generality, the norm | • | is sub-multiplicative. 
Since A is bounded, there exists M > 0 such that, for every A G A, j G |1, IV], and 
t G R, we have |Aj(t)| < M. Let A G A. For every n G N^, let be a set of 

measure zero such that < f{L-n) holds for every t G R \ N[n]- Let N = 

UnGN" which is of measure zero. If n G and t G {L ■ n — Lmax, L ■ n)\N, 
then 


0 


L,A,A 

[n],t 


< 


IMJ 

L-n—Li<.t 


^L,A,A 

“[n-ej],* 


Aj] (t — L ■ n -G Lj) 


< NM 


■j -i*- 

E 

li]eJ 


f{L-n- Lj) < C 


max 


/(s), 


where C = N'^M. 


□ 
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As an immediate consequence of the previous proposition and Theorem 3.21, 
we have the following theorem, which improves Theorem 3.21 by replacing (0, A)- 
exponential type by (S, A)-exponential type. 


Theorem 3.25. Let A G (Klj_)^ and A be a bounded shift-invariant subset of 
Ald(C)^). For every L G V+(A), Sg(L,A) is of CE., A)-exponential type 7 if 
and only if it is of {Q, K)-exponential type 7. 

As a eonsequence, for every L G T+(A), if Yis{L, A) is of {E, A)-exponential type 
7 then, for every p G [l,+oo], it is of exponential type 7 in Conversely, for 
every L G W+{A), if there exists p G [l,+oo] sueh that 'Fs{L,A) is of exponential 
type 7 in X^, then it is of iE, A)-exponential type 7. Finally, for every L G W+(A) 
and p G [1, +00], 


Xp{L, A) 


lim sup sup ess sup 

Z/'ii— 


In 

Si.AA 


[n]7 


L • n 


(24) 


3.3.3. Arbitrary switehing. We consider in this section A of the type A = L°°(R, *8) 
with ® a nonempty bounded subset of Aid{C)^. In this case, F,s{L,A) corresponds 
to a switched system under arbitrary 18-valued switching signals (for a general 
discussion on switched systems and their stability, see e.g. [20,28] and references 
therein). 

Motivated by formula (12) for j, we define below a new measure of the 
asymptotic behavior of Y,s{L,A). For this, we introduce, for A G (Rl^)^ and 
X G R+, 

£(A) = {A • n I n G N^} and £a;(A) = £(A) fi [0, x). (25) 


Definition 3.26. We define 


_1 

X 


p(A, 18) = limsup sup 

X—>-+oo 

xGC{A) for r^Cx{A) 


|n|l 




A-p„(fe) 

Vk 




Note that /r(A, 18) is independent of the choice of the norm j-j and /r(A, 18) = 
p{A, 18). The main result of this section is the following. 

Theorem 3.27. Let A G (R^)'^, L G V+(A), ^ be a nonempty bounded subset of 
Md{C)^, A = L°“(R, 18), and p G [l,-|-oo]. Set mi = minjg|i_ 7 v| ^ and m 2 = 
maxjg|i jv] ^ *//7(A, 18) < 1 , and mi = maxjg|i jv] ^ and m 2 = minjg|i_jv] ^ if 
/r(A, 18) > 1. Then the following assertions hold: 

i. Ap(L, .4) < mi ln/i(A, 18); 

ii. if L G bF+(A), then m 2 Xp{A,A) < Xp{L,A) < miXp{A,A); 

iii. Ap(A,Al) = ln/r(A, 18). 


Proof. Notice that (ii) follows from (i) and (iii) by exchanging the role of L and A, 
since A G V+(L) for every L G 14 +(A). 

Let us prove (i). Since min^gii ^v] ^ E E niciXjg|i iv] ^ for every n G 
\ {0}, it suffices to show that, for every e > 0, there exists C > 0 such that, for 
every A G Al, n G \ {0}, and < G R, we have 
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By definition of /i(A, S), there exists Xq € /i(A) such that, for every x S C[K) 
with X > Xq, we have 


sup 

B^e<3 

for r£Cx{A.) 


A-n=a: 


<(M(A,»)+er. 


Since fB is bounded, the quantity 


C' = max sup 
x^CxqX) B^e<s 

for r^Cx{A) 


E E 


^’^Vn k—1 

An—x 


is finite. Setting C = max{l, C", C"(^(A, 55)+e) we have, for every x G B(A), 


sup 

B^e<s 

for rGCxiA) 


E E n s™"-'*’ 

A-n—r 


< C + e)\ 


(27) 


Define (pL ■ B(A) —>• C{L) by (Pl{A ■ n) = L ■ n. This is a well-defined function 
since L G V+(A). Let A G n G \ {0}, and t G M. By Proposition 3.12, 


"L,A,A 

“[n],t 


Yl AvAt - L ■ Pv{k)). 

n'G[n]nN''JJGVn/ fc=l 


(28) 


For r G £A.n(A), we set i?'’ = A{t — pl{t)) G 55. Thus, for every n' G [n] fl N^, 
V G Vn', and k G |1, |n'|^], we have, by definition of pL, 

^A.p„(fe) ^ -pl{A- Pv{k))) =Ay^{t-L- Pv{k)) . (29) 


We thus obtain (26) by combining (27), (28) and (29). 

In order to prove (iii), we are left to show the inequality ln/i(A, IB) < Ap(A,A). 
Let X G £(A) and A° G IB. For r G Cx{X), let S’" G IB. We define 


c 


— min 1 2/1 

2 yi,V2^Cx{A) 
yi/y2 


2 / 2 I > 0. 


Let A = [Ai 


.., Ajv) G A be dehned for t G K by 

{ , if m G is such that A • m < x 

^Am 

and < G (-A • m - C,-A • m + C), 
A°, otherwise. 


The function A is well-defined since the sets (—A • m — — A • m -|- C) are disjoint 
for m G with A • m < a:. For every n G with A • n = X, every v G Vn, 
t G (-CO, and k G |1, |ny, we have 

Ax, (t-A.p„(fc)) = <'’’'^"\ 

and then, for every n' G with A • n' = a;, we have 




k—1 

A-n—x 


l“li 

E E n - A • Pv{k)) 

nGN-''' ’'GVn k=l 
A-n—x 


GA,A,A 
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Hence, for every n' € with A • n' = X, we have 


veVn k=l 

A-n—a: 


< 


sup ess sup 
AeA tGR 



Since this holds for every choice of S’" G *B, r G and x G >C(A), we deduce 

from (24) that ln/i(A, IB) < Ap(A,^). □ 


Remark 3.28. Since = ^(A, 55), one has Ap(A,Al) = Ap(A, 55)). 

As regards exponential stability of Tis{L,A), we deduce from the previous theo¬ 
rem and Remark 3.16 the following corollary. 

Corollary 3.29. Let A G (K.^)^, ^ be a nonempty bounded subset of Md(C)^, 
and A = L°°(K., IB). The following statements are equivalent: 

_i. MA,») < 1; 

ii. E 5 (A,AI) is exponentially stable in for some p G [l,-|-oo]; 

ill. T,s{L,A) is exponentially stable in for every L G V+(A) and p G [1, -l-oo]. 

Moreover, for every p G [1, -l-oo], 

Ap(A, A) = inf{i/ G M I p(A, 55_,y) < 1}, 
where 55_, = {{e-'^^Bi ,..., | (Ri, • ■ •, Riv) e $}• 

Corollary 3.29 is reminiscent of the well-known characterization of stability in the 
autonomous case proved by Hale and Silkowski when A has rationally independent 
components (see [4, Theorem 5.2]) and in a more general setting by Michiels et al. 
in [21]. In such a characterization, (!,...,!) is assumed to be in H(A) and /r(A, IB) 
is replaced in the statement of Corollary 3.29 by 

Phs(A, A) = max P 
(ei....,e„)GU(A) 

where H(A) is the image of H(A) by the canonical projection from onto the 
torus (R/27rZ)^. (Notice that H(A) is compact since the matrix B characterizing 
H(A) in Proposition 3.9 has integer coefficients.) 

We propose below a generalization of pHs(A,A) to the non-autonomous case 
defined as follows. 

Definition 3.30. For A G (R^)^, IB a nonempty bounded subset of , and 

£(A) given by (25), we set 

X 
n 

E • 

Let us check in the next proposition that pns actually extends pns ■ 

Proposition 3.31. Let A = (Al,..., An) G A4d{C.)^ and 55 = {A}. Then one 
has phs(A, IB) = phs(A, A). 


Mhs(A, IB) = limsup sup 


n—>-+oo 


sup 


(Si 
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Proof. One has 


max _ p 

(0i,...,e;v)6R(A) 



max _ lim 


lim sup 


[P-'j 

iP-j 


”^+°°(ei,....e«)GY(A) „g|i,Ar|nfe=i 


where the second equality is obtained as consequence of the uniformity of the 
Gelfand limit on bounded subsets of A4d(C) (see, for instance, [13, Proposition 
3.3.5]). □ 


In the sequel, we relate /rHs(A,f8) to a modified version of the expression (24) 
of Xp(L, A). 


Definition 3.32. For L G 1 + (A) and A a set of functions A : K. —>■ Ald(C)'^, 
define 


Ahs(A,.4) 


lim sup sup ess sup 
|nq->-+ooAGA tGR 


In 

"L,A,A 


[n],t 


n 


we 


Remark 3.33. Since Lminjnj^ < L ■ n < Lmaxlnj^ for every L G V+(A) and 
n G N^, it follows immediately from (24) that, for every p G [1, +oo], 

^min Ap {L,A) < Ahs(A,^) < ^maxAp {L,A), if Ap(L,Al) >0, 

^maxAp (L,A) < AHs(.b,Al) ^ ^min Ap(L,Al), if Ap(L,Al) <0. 


In particular, the signs of Xhs{L,A) and Xp{L,A) being equal, they both charac¬ 
terize the exponential stability of '£s{L,A). 


^)^, ^ be a nonempty bounded subset of X4d{C)^, 
b+li 


z G 


Z(A)\{0}} *//rHs(A,») 


< 1 


Theorem 3.34. Let A G 

and A — L°°(R, 25). Set m = inf|l,p—p 

and m = sup|l,pA-p zG .Z(A)\{0}| if /iHs(A, 25) > 1. Then the following 
assertions hold: 

i. for every L G 14(A), Ahs(A,AI) < TOln/iHs(A, ®); 
ii. «/ (1,..., 1) G P(A) and L G IF4(A), one has AHs(.b,4) = ln/iHs(A, 25). 


Proof. We start by proving (i). It is enough to show that, for every £ > 0 small 
enough, there exists C > 0 such that, for every 4 G 4, n G \ {0}, and t G R, 
we have 

4n]y I < ^(1 + |n|i) (mhs(A, ») + . 

Let L G 14 (A) and £ > 0 be such that phs(A, IB) -I- £ < 1 if /rHs(A, IB) < 1. We 
can proceed as in the proof of Theorem 3.27 to obtain a finite constant Co > 0 such 
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that, for every n gW, 


sup 


sup 


for reiZnAmaxC^) 


j;G[l,Afl" k=l 


^■Pv{k) 


< Co (mhs (A, ») + £)”. 

(30) 

Let A G A, t G M, and (fL be as in the proof of Theorem 3.27. For r G /iraA„ax(A), 
we set = A{t — ifL{r)), and similarly to the proof of Theorem 3.27, (29) holds 
for every v G lljiV]" and k G |l,n|. Thus (30) implies that, for every n G N* and 
0 gL(A), 


k=l 


Since 


< Co (^hs(A, 25) + £)" . 


|n|i 


Yi^vAt-L-Pv{k))e^^-'‘ = X! 

DG[l.Af]" k=l 


nGN™^ veVn fe=l 

|n|i=rt 


|n|l 


= E n^-ut-A-p.(fc)) 


nGN^ i;GVnfe=l 

|n|i=rt 


= E 


^in-9'^L,A 
" “n.i ? 


|n|i=n 


we obtain that, for every n G N* and 0 G F(A), 


E 

nGN" 

|n|i=ra 


^in-d'^L.A 

“n.i 


< Co (/xhs(A, 18) + e)" . 


(31) 


Following Proposition 3.9, fix h G |l,iV] and B G AiN,h(A) with rk(i?) = h 
such that A = Bio for io G (RlJ.)^ with rationally independent components. Let 
M G A1a(R) be an invertible matrix such that £q = Mei, where ei is the first 
vector of the canonical basis of K^, in such a way that A = BMei. For n G N, we 
dehne the function /„ : —>• Md{A by 


fniv) = E 


^in.- B M y '^L. A 
“n,i ■ 


|n|l=" 


We claim that, for every no G N^, 

1 f 


lim 


fi -!-+00 { 2R )^ 


fnAy 


y: 


(32) 


nG[no]nN' 

|n|i=n 


Indeed, we have 
1 f 


(27?)'* J[-R^R]h 


fn{vy 


-ixiQ-BM. 


•'diy= Y. 


"L,A 


J{n-no)-BMu 


nGN" 

|n|i=r: 


(27?)'^ Jl-R,R]>^ 


dv. 
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If n S is such that A • n = A • no, then A • (n — no) = 0, and therefore 
n — no £ C y(A)-*- = (Rani?)-*-. One gets (n — no) • BMv = 0 for every 

1 / £ implying that 


1 

{2RY 




dv = 1. 


[-R,R]'' 


If now A • n 7 ^ A • no, set = A • (n — no), which is nonzero. Then 


I 

{2RY 


Jin-no)- 


[-R,R]^ 




sin(^i?) 


-^0, 

R —>^ + CX5 


which gives (32). 

We can now combine (31) and (32) to obtain that, for every n £ N* and no £ 

N^\{0}, 


E 

nG [noInN”^ 

|n|i=n 


^L,A 

“n.t 


< C'o(/^Hs(^i ®) + sY' 


Set mo = sup|l,||ip zG Z(A)\{0}| and notice that, since Z{K) = —Z{A), 
one has = inf |l, z £ Z{K) \ {0}|. We claim that, if n, no £ and 

A • n = A • no, then ^ |no|^ < |n|^ < mo |no|^. Indeed, let z = n — no £ Z{K) and 
ni = no — z_ £ N^. Then one has 


1 

-, mo 

mo 

Hence, for every no £ \ {0}, 


Hi ^ l^+li + l^ili ^ 
Inoli |z-L+|niL 


^L,A,A 

“[no],t 


+ 00 

E E -n:" 


"■=OnG(no]nN" 

|n|i=n 


nG 


E 


E 


ne[no]nN^ 

^ ^ |n|j^=n 


^L,A 


and we conclude that 

E Co(ms(A,S)+e)’^<C(I+|no|i)(MHs(A,$)+e)'"l"“ls 

"*o|nolil 


for some C > 0. This concludes the proof of (i). 

Suppose now that (I,..., 1) £ R(A). Then |z+|^ = \z-\i for every z £ Z{A), and 
hence (i) yields Xhs{L,A) < ln^Hs(A, IB) for every L £ V+(A). We claim that it is 
enough to prove (ii) only for L = A. Indeed, assume that Ahs(A,. 4) = ln^Hs(A, fB). 
In particular, 

•^hs(A,AI) < Ahs(A,AI) (33) 

for every L £ V+(A). Since A £ V+{L) if L £ 1T+(A), by exchanging the role of L 
and A in (33), we deduce that Ahs(A,AI) = Ahs(A,AI) for every L £ W+(A), and 
hence (ii). 

Let n £ N* and B’' £ 05 for r £ CnA^^^iA). As in the argument for (iii) in 
Theorem 3.27, there exist C > 0 and a function A : R —>■ A4d{C)^ such that, for 
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every v S |1,-/V|", t S (—CiO) and k G |l,n|, we have 

n 

A,, {t-A- p.(fc)) = and ^ J] ^ 

t;G[l,AIl" fe=l nGN" 

|n|i=rt 


Denote Z+ = {[n] G Z | [n] fl ^ 0}. Since (1,..., 1) G D(A), one deduces 
that, if n, n' G are such that n si n', then e*"'® = e*" for every 9 G D(A) and 
|n|^ = |n'|^. We set |[n]|^ = |n|^ for every n G N^. Then 


E 




nGN" 

|n|i=n 


E 


E 


in •d'^'A,A 

^ “n',t 


[n]^Z+ n'G[n]nN^ 
|[n]|i=r!, 


E 


m-0^A,A,A 


[n]G2+ 

|[n]|i=n 


We clearly have #{[n] G Z+ | |[n]|^ = n} < #{n G | |n|^ = n} = < 

(n + 1)^“^, and we get that, for every 9 G V{A) and n G with |n|^ = n, 



n 

~ 


E 

JJ- 5A.p„(/c)gze„, 

N-l 

< (n -1-1) " sup ess sup 

SA,A,A 

fsiiwr 

k^l 

AeA teR 



Since the above inequality holds for every choice of S’' G fB, r G B„a^^^(A), n G N*, 
we deduce that ln/rHs(A, ®) < Ahs(A,.4). This concludes the proof of Theorem 
3.34. □ 


The next corollary, which follows directly from the above theorem and Remarks 
3.16 and 3.33, generalizes the stability criterion in [4,21] to the nonautonomous case 
(see Proposition 3.31). 

Corollary 3.35. Let A G ^ be a nonempty bounded subset of AidiC)^, 

and A = L°°(R, *8). Consider the following statements: 
i- /iHs(A, *8) < 1; 

ii. E 5 (A,AI) is exponentially stable in for some p G [l,+oo]; 

ill. IisiL,A) is exponentially stable in for every L G V+(A) and p G [1, +oo]. 

Then (i) => (Hi) => (ii). If moreover {1,... A) ^ V(A), we also have (ii) 

(i) and, for every p G [1, +oo], 

Ap(A, A) = inf{u G K | AtHs(A, *8-,.) < 1}, 
where = {{e-'^^Bi,..., e-'^^ B^) |(Si,..., Sjv) e »}• 

4. Transport system. For L = (Li,... ,Ln) G and M = (wy )ijg|i_Ar| : 

R —^ A4Ar(C), we consider the system of transport equations 

{ du' du' 

x) + -^{t,x) = 0, ie |1, A^l, < G [0, +oo), a; G [0, L,], 

N ^ 

Ui{t,0) = '^mij{t)uj{t,Lj), i G [1,A^1, t G [0,+oo), 
i=i 

(34) 

where, for i G |1, N}, Ui{-, •) takes values in C. 

The time-varying matrix M represents transmission conditions and in particular 
it may encode an underlying network for (34), where the graph structure is deter¬ 
mined by the non-zero coefficients of M. When no regularity assumptions are made 
on the function M, we may not have solutions for (34) in the classical sense in 
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CHK+ X [0,L*]) nor in C'°(R+, Li], C)) n C'i(IR+, LP([0, L*], C)). We thus 

consider the following weaker definition of solution. 

Definition 4.1. Let M : K —>■ A4n{C) and Ui^o : [0,Li] —>■ C for j S |1,7V|. We 
say that (■Ui)iG|i,A] is a solution of Tir{L,M) with initial condition (ni,o)iG[i,Afl if 
Ui : M+ X [0, Li] —>• C, i G |1, NJ, satisfy the second equation of (34), and, for every 
^ 0, X G [0, Li], s G [— min(x, t), Li—x], one has X“t“s) — Ui(t^ x) 

and Ui(0,x) = Uifi{x). 

4.1. Equivalent difTerence equation. For i G |l,iV| and M : K — 

define the orthogonal projection Pi = aej and set ^i(-) = M{-)Pi. Consider the 
system of difference equations 

N 

v{t) = '^Aj{t)v{t - Lj). (35) 

i=i 

This system is equivalent to (34) in the following sense. 

Proposition 4.2. Suppose that (ui)iG|i,A] is a solution of (34) with initial condi¬ 
tion (■Ui,o)iG[i,A'l '■ [~^max, +oo) —>■ C'^ be given for i G |1, iV] by 

{ 0, f G [ Lniax, Li) , 

Mi,o(-t), ift€[-Li,0), (36) 

Ui(t,0), ift>0. 

Then v is a solution of (35). 

Conversely, suppose that v : [—Lmax, +oo) —>■ is a solution of (35) and let 

(Mi)ig|i^Ar] given for i G |1, A^jj, t > 0 and x G [0,Li] by Ui{t,x) = Vi{t — x). 
Then ('ni)iG|i,Ar] is a solution of (34). 

Proof. Let (wijigii.Ari be a solution of (34) with initial condition ('ai,o)iG[i.Afl and 
let V : [—Lmax, +oo) —>■ be given by (36). Then, for t > 0, 

N 

Vi{t) =Ui{t,0) = 'y^mij[t)uj{t,Lj), 

1=1 

and, by Definition 4.1 , Uj(t, Lj) = Vj(t — Lj ) since Uj(t, Lj ) = Uj(t — Lj,0) iit > Lj 
and Uj(t,Lj) = Ujp{Lj — t) if 0 < t < Lj. Hence Vi(t) = J2j=i ~ ^l) and 

thus v{t) = J2f^i Aj{t)v{t - Lj). 

Conversely, suppose that v : [—Lmax, +oo) —>■ is a solution of (35) with 

initial condition vq and let (Mi)iG|i,A] be given for i G |l,iV], t > 0 and x G [0, Li] 
by Ui{t,x) = Vi(t — x). It is then clear that Mi(t + s,x + s) = Ui{t,x) for s G 
[— min(x, t), Li — x], and, since Vi(t) = ~ -bj), 

N N 

Ui{t,0) = Vi{t) = '^mij{t)vj{t - Lj) = y^mij{t)uj{t,Lj), 

1=1 1=1 

and so (ui)ig|i_jv] is a solution of (34). □ 

The following result follows immediately from Proposition 3.2. 

Proposition 4.3. Let Ui^ : [0, Li] ^ <C for i G fl, NJ and M :'R ^ Mn{<C) . Then 
Et-(L, M) admits a unique solution (■Ui)iG|i,A ]7 Ui : M+ x [0,Li] —i C fori G |1, Ai], 
with initial condition (iti,o)iG[i.AiT 


STABILITY OF DIFFERENCE EQUATIONS AND APPLICATIONS 


25 


4.2. Invariant subspaces. For p G [1, +oo], consider (34) in the Banach space 


N 




endowed with the norm 


l|M|lp = 


N 

El 


i=l 


^i|lLP([0,Li].C) 


1/p 




if p G [1, +oo), 
ii p = + 00 . 


It follows from Proposition 4.2 and Remark 3.4 that, if M G AIjv(C)) and 

uq G Xp, then the solution t i—^ u{t) of T,t{L,M) with initial condition uq takes 
values in XJJ for every t > 0. 

In Section 5, we study wave propagation on networks using transport equations 
via the d’Alembert decomposition. For that purpose, we need to study transport 
equations in the range of the d’Alembert decomposition operator, which happens 
to take the following form (see Proposition 5.3). For r G N and R G A4r,N{C) with 
coefficients pij, i £ |I,r], j G |I,fV|, let 


Yp(i?) = < M = (mi, ..., un) G Xp 


N 


Y* e |I,r], Vpij / Uj{x)dx = 0 > . 


This is a closed subspace of Xl), which is thus itself a Banach space. 


Remark 4.4. Let r G N, A G A^r.Ar(C), and M G A4iv(C)). Note that, if 

1 < P < < + 00 , yq{R) is a dense subset of Yp(i?) since X^ is a dense subset of Xp. 

As a consequence, by a density argument, Propositions 3.14 and 4.2, one obtains 
that, if Yp(i?) is invariant under the flow of Sp(L,M) for some p G [l,+oo], then 
Yq(i?) is invariant for every q G [1, +oo]. 


The following proposition provides a necessary and sufficient condition for Yp(ii) 
to be invariant under the flow of (34). 


Proposition 4.5. Let r G N, i? G AlrAr(C), (uio)iGliAfl ^ Yp(i?), and M G 
A4jv(C)). Then the solution u = (ni)ig|i_Ar| of T,riL, M) with initial condi¬ 
tion (ui_o)iGli,Afl belongs to Yp(i?) for every t > 0 if and only if 

R{M(t) — ldN)w(t) = 0 

for almost every t>0, where w = (rei)jg|iyv] eind Wi(f) = Ui(t,Li). 


Proof Let v : [— Lmax, +oo) — >■ be the solution of (35) corresponding to u, 

given by (36), and let w = (wi)iG[i,Afl be defined by Wi{t) = Viit — Li) = Ui(t,Li). 
Let A = (Ai)jg|i_p| be given for i G |I,r] by Xi{t) = Pd fo' Uj{t,x)dx. Since 

Ai(0) = 0, we have 


N 

(^) = E! Pd 

N 

= E 

j=i 


r^i 

r^j 


/ Uj{t,x)dx — 

Ujfl{x)dx 


Jo J 

0 


r^3 

r^j 

- 

1 

H 

H 

1 

/ '>Jj{-x)d 

X 

/o 

Jo 
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N 

= X! 

i=i 


j{s)ds— / Vj{s — Lj)da 


t — Lj 


= ^P^3 / iv3is)-Vj{s-Lj))ds 

3 = 1 -^0 

N t ( N \ 

= '^P ^3 y] rnjkis)vkis - Lk) - Vj{s - Lj) ds 

J = 1 -^0 \fc=l / 


N .t N 


"^'^Pi 3 / y (mjkjs) - 5jk) Vkjs - Lk)ds, 

j = l “'0 fc=l 

so that \{t) = fg R(M(s) — ld]v)w(s)ds. The conclusion of the proposition follows 
immediately. □ 


Definition 4.6. Let L G (R!j_)'^ and At be a subset of At iv(C)). We denote 

by Inv(Af) the set 

Inv(Af) = {i? G Afr,Ar(C) | r G N, Yp(i?) is invariant under 
the flow of T,r{L, M), VM G Af, Vp G [1, +oo]}. 

4.3. Stability of solutions on invariant subspaces. We next provide a defini¬ 
tion for exponential stability of (34). 

Definition 4.7. Let p G [l,+oo], L G Af be a uniformly locally bounded 

subset of L^^(M.,M.n(C)), and R G Inv(A4). Let ST-(L,Af) denote the family of 
systems 'E,t{L,M) for M G A4. We say that St-(L,A 1) is of exponential type 7 in 
Yp(i?) if, for every e > 0, there exists A > 0 such that, for every M G M and 
uq G Yp(i?), the corresponding solution u of M) satisfies, for every t > 0, 

We say that (L, A4) is exponentially stable in Yp(i?) if it is of negative exponential 
type. 

The next corollaries translate Propositions 3.19 and 3.20 into the framework of 
transport equations. 

Corollary 4.8. Let A G L G V1)_(A), and A4 be a uniformly locally bounded 

subset of A4 n(C)). Suppose that there exists a continuous function /:!&—>■ 

K!j_ such that, for every M G M, n G N^, and almost every t G (L-n—Lniax,L-n), 
(17) holds with A = (Ai, ..., Av) given by Ai = MPi. Then there exists a constant 
C > 0 such that, for every M G Af, p G [l,+oo], and uq G Xp, the corresponding 
solution u of Yir{L,M) satisfies 

||w(t)||p < (^(t + 1)^”^ max /(s)||Mo||p, Vt > 0. 

sG[t ■L'm.a.xi*] 

Proof. Let C > 0 be as in the Proposition 3.19. Let M G Af, p G [l,+oo], 
wq G Xp, and u be the solution of T,r{L,M) with initial condition uq. Let v be 
the corresponding solution of (35), given by (36), with initial condition vq. Notice 
that lluollp = ||uo||p and, for every t > 0, ||u(t)||p < ||ut||p. By Proposition 3.19, we 
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have, for every t > 0 , 

II^WIIp < Iktilp < C'(f + 1 )^-^ max /(s)||po||p 
= C(t + l)^-^ max /(s)||uo||p, 

which is the desired result. □ 


Corollary 4.9. Let A € L € W+(A), Ai be a uniformly loeally bounded 

subset of AfAr(C)), and / : R —?> be a continuous function. Sup¬ 

pose that there exist M G A4, Uq G N^, and a set of positive measure S C 
(L ■ no — Amaxj dy ■ no) such that, for every t G S, (19) is satisfied with A — {Ai ,..., 
An) given by Ai = MPi. Then there exist a constant C > 0 independent of 
f, an initial condition uq G X^, and t > 0 such that, for every p G [l,+oo] 
and R G Inv(Ad), the solution u of Y,r{L, M) with initial condition ug satisfies 
u{s) G Yp(i?) for every s > 0 and 

\\u{t)\\p>C min /(s)||uo||p. 


Proof. As in Proposition 3.20, since L G W-|_(A), we can assume for the rest of the 
argument that A — L. 

Let C > 0 be as in Proposition 3.20. We construct an initial condition vq G Xp 
as follows: choose to and jo as in Proposition 3.20 and verifying in addition to ^ 
L • no — Ljg. Then pick 5 > 0 as in Proposition 3.20 and satisfying in addition 
i5 < |to — A ■ no + Ljg \ and 6 < Linin/2. Next, take fj, G A°°(R,M) as in Proposition 
3.20 and satisfying in addition fi(s)ds = 0. Finally, consider the initial condition 
z;o(s) = p(s — to -G L ■ no)ejg. As in (22), we still obtain that the solution v of (35) 
with initial condition vg satisfies, for p G [1, +oo], 

Ikto+illp > lktollLP([-5,5],c'V) > C . , min fis)\\vo\\ . (37) 

^ ' sG[to + l5-Lmax,to + o] 


Let u be the solution of (34) corresponding to v, in the sense of Proposition 4.2, 
and denote by ug = (wLo)iG[i,Ail its initial condition. Since Ui^o{x) = Vi{—x), we 
have Ug G Ai],C). Furthermore, Ui^o = 0 for f ^ jg and Ujgfi{x) = 

Vjg (—x) = fi{L ■ Uq — tg — x). By definition of S, we must have either {L ■ ng — tg — 
S, L-ng — tg-GS) C [0,Ljo] or {L-ng — tg — d,L-ng — tg-G6)r\[0,Ljg] = 0, but the latter 
case is impossible since we would then have Ujg^ = 0 , and thus u(s) = 0 for every 
s > — Aniax, which contradicts (37). Hence {L ■ ng — tg — 6, L ■ ng — tg -G S) C [0, Ljg] 
and 

rLjg rS 

/ Ujgfi{x)dx = / p,{x)dx = 0 . 

Jo ’ J-s 

We thus have clearly ug G Yac,(R), and in particular u(s) G Yp{R) for every s > 0 
and p G [1, +oo]. Furthermore, Hiiollp = ll^oUp and, for p G [1, +oo). 


j-s j-s N 

\\'^to\\lp([- 5 ,s],c^)= + s)|pds = / ^ |u,(to+ s,0)|^ds 

J-s J-s 


1-5 

I-2S JV N Li 

^ |ui(to+ ^, s)|^(is < ^ / |ui(fo + <j, s)|^ ds 




i=l 


i=l ■ 


= IKio + d)||^ 
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with a similar estimate for p = +c». Hence, it follows from (37) that, for every 
p e [l,+oo], 

\Wit)\\p>C min /(s)||uo||p 

sG [t ^max 

with t = to + S. □ 

As a consequence of the previous analysis, we have the following result. 

Theorem 4.10. Let Ai he a uniformly locally hounded subset of L-(R,Af^(C)), 
A e (K;)^, and A = {A = {Ai,...,An) : M ^ Mjv(C)^ | Ai = MPi,M G 
M.}. For every L G V+(A), if Yis{L,A) is of {Q, A)-exponential type 7 then, for 
every p G [l,+oo] and R G Inv(Af), T,r(L,Ai) is of exponential type 7 in Yp{R). 
Conversely, for every L G W+{A), if there exist p G [l,+oo] and R G Inv(Al) 
such that AI) is of exponential type 7 in Yp(R), then Yis{L,A) is of (0, A)- 

exponential type 7 . 

It follows from Theorem 4.10 that the exponential type 7 for Tjr{L,Ai) in Yp{R) 
is independent of p G [1, + 00 ] and R G Inv(AI). When Ai is shift-invariant, thanks 
to Theorem 3.25, one can replace (0, A)-exponential type by (S, A)-exponential 
type for '£is{L,A) in Theorem 4.10. 

Assume now that AI = where IB is a bounded subset of AlAr(C). 

Let A = {A = (Ai,...,AAr) : R ^ A4Ar(C)^ | A, = MP„ M G Ai}, i.e., A = 
L“(K,2t) where 21 = {A = {Ai,...,An) G A4jv(C)^ | A, = MPi, M G $}. We 
can thus transpose the results from Section 3.3.3, and in particular Corollary 3.29, 
to the transport framework. 

Corollary 4.11. Let A G (R^)-^, ^ be a nonempty hounded subset 0 /A4Ar(C), 
AI = L°°(R, 05). The following statements are equivalent. 

i. E^(A,A4) is exponentially stable in Yp(i?) for some p G [l,-|-oo] and R G 
Inv(Al). 

ii. 'Fr{L,Ai) is exponentially stable in Yp(i?) for every L G V+(A), p G [l,-|-oo], 
and R G Inv(A4). 

Remark 4.12. In accordance with Remark 3.28, the exponential stability of the 
system Sr(A, A4) is equivalent to that of Sr(A, L°“(R, IB)). 

5. Wave propagation on networks. We consider here the problem of wave prop¬ 
agation on a finite network of elastic strings. The notations we use here come 
from [ 11 ]. 

A graph C/ is a pair (V,f), where V is a set, whose elements are called vertices, 
and 

E C {{q,p} I q,p G V, q^p}. 

The elements of £ are called edges, and, for e = {q,p} G £, the vertices q,p are 
called the endpoints of £. An orientation on Q is defined by two maps a,ui : £ ^ V 
such that, for every e € £, e = {a(e),a;(e)}. Given q,p G V, a path from q to p is 
a n-tuple {q = qi,... ,qn = p) G V" where, for every j G |1,n — 1], {qj, qj+i} G £. 
The positive integer n is called the length of the path. A path of length n m. Q 
is said to be closed if qi = simple if all the edges {qj,qj+i}, j G |l,n — 1|, 
are different; and elementary if the vertices gi,..., are pairwise different, except 
possibly for the pair (gi,g„). An elementary closed path is called a cycle. A graph 
which does not admit cycles is called a tree. We say that a graph Q is connected 
if, for every g,p G V, there exists a path from q to p. We say that Q is finite if V 
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is a finite set. For every q € V, we denote by £q the set of edges for which q is an 
endpoint, that is, 

£q = {e € £ \ q € e}. 

The cardinality of £q is denoted by Uq. We say that g G V is exterior if £q contains 
at most one element and interior otherwise. We denote by Vext and Vint the set 
of exterior and interior vertices, respectively. We suppose in the sequel that Vext 
contains at least two elements, and we fix a nonempty subset Vd of Vext such that 
Vu = Vext\Vd is nonempty. The vertices of Vd are said to be damped^ and the vertices 
of Vu are said to be undamped. Note that V is the disjoint union V = Vint U Vu U Vd. 

A network is a pair (C/, L) where G — (V, £) is an oriented graph and L — {Le)eeS 
is a vector of positive real numbers, where each is called the length of the edge 
e. We say that a network is finite (respectively, connected) if its underlying graph 
G is finite (respectively, connected). If e € £ and u : [0, Lg] —>• C is a function, we 
write u(a(e)) = u(0) and u(uj(e)) = u(Le). For every elementary path (gi,... ,gn), 
we define its signature sif—J-l— 1 , 0 , 1 } by 

{ 1 , if e = {gi,gi+i} for some i G |l,n — I] and a(e) = qi, 

— 1, if e = {qi,qi+i} for some * G |l,n — I] and a(e) = g^+i, 

0 , otherwise. 


The normal derivatives of u at a(e) and uj(e) are defined by -^{o.{e)) = —fy(0) 
and|^(..(e)) = ^(Le). 

In what follows, we consider only finite connected networks. In order to sim¬ 
plify the notations, we identify £ with the finite set |1,A^], where N = ff£. We 
model wave propagation along the edges of a finite connected network {G, L) by N 
displacement functions Uj : [0,-|-oo) x [0,Lj] —>■ C, j G |1, A^], satisfying 


T,u:{G,L,f]) : 



Uj{t,q) 


d^Uj , 

Uk{t,q), 


0 , 


j G |1, NJ, t G [0, -boo), x G [0, Lj], 
q G V, j,k G £q, t G [0, -boo), 

g G Vint) t G [0, -boo), 


(38) 


UU ■ 

Q e Vd, j G£q, tG [0,+00), 
0, g G Vu, j G £q, t G [0, -boo). 


Each function rjq is assumed to be nonnegative and determines the damping at the 
vertex g G Vd- We denote by rj the vector-valued function rj = {r]q)q^Vd- 

Remark 5.1. Let {G,L) be a finite connected network with £ = |1,A^] and 
(ai,tui), (a 2 ,W 2 ) be two orientations of G- If (^^j)jG[i,Ail satisfies (38) with ori¬ 
entation (ai,a;i) and (wjOjgIi.ivI is given by Vj = Uj if ai(j) = a2(j) and Vj{x) = 
Uj{Lj — x) otherwise, we can easily verify that ('yj)je[i.Ail satisfies (38) with ori¬ 
entation (a 2 ,W 2 ). Hence the dynamical properties of (38) do not depend on the 
orientation of G- 
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For p S [1, +oo], consider the Banach spaces LP{Q, L) = O^i -^^([0; Lj\,C) and 


N 


L) = I (ui,..., Wat) G J] 

Uj{q) = Ukiq), Vg G V, Vj, fc G £q\Uj{q) = 0, Vq G Vu, Vj G £q} , 


endowed with the usual norms 


N 


\Lp(g,L) 


E 11^*11 


LP([0,Li],C) I ’ 


^ 2=1 




if p G [1, +oo). 
if p = +00, 


N 


l^llwo'P(e,L) 


E ll^'illinio.i.i.c) ’ if P e [1, +oo), 


^ 2=1 


> if p = +00- 


^eli.ATl 




We will omit {G,L) from the notations when it is clear from the context. 

Let Xp = Wq’^ X LP, endowed with the norm Ij-Hp defined by 

+ MlnG,L ))" ’ if p e [1, +oo), 
max , if p = +oo, 

and, for every t G M, define the operator A{t) by 

D{A{t)) = |(zr,u) G n W2'P([0,L,],C)j x 

Vjiq) = Vg G Vd, Vj G £q; E ^(9) = 0, Vg G Vi„t ^ , 


„ drij 
jee, ^ 


A{t) 


(39) 


V J \u 

One can then write (38) as an evolution equation in as 

U{t) = A{t)U{t) 

where U = (w, fy)- 

5.1. Equivalence with a system of transport equations. In order to make a 
connection with transport systems, we consider, for p G [1, +oo], the Banach space 

2N 

3 = 1 

where for j G |l,fV|. 


Definition 5.2 (D’Alembert decomposition operator). Let T : X^ —>■ XJ be the 
operator given by T(m, v) = /, where, for j G |1, iV|, cc G [0, Lj], 

/ 2 j _ 1 (x) = u' (Lj - x) + Vj (Lj-x), f 2 j{x) = u'^{x) - Vj{x). (40) 
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In order to describe the range of T, we introduce the following notations. Let 
r £ N be the number of elementary paths (gi,..., q„) in Q with qi = qn or qi,qn £ 
Vu- The set of such paths will be indexed by |1, r|. We denote by Si the signature 
of the path corresponding to the index i £ |1, r|. We define R £ AIr, 2 iv(C) by its 
coefficients pij given by 

= Pz, 2 j = s^{j) for i £ |l,r], j £ [ 1 , 

We then have the following proposition. 

Proposition 5.3. The operator T is a bijection from Xf to Yp{R). Moreover, T 
and T~^ are continuous. 


Proof. Let {u, v) £ and let / = T{u, v) £ XJ. Let (qi,..., qn) be an elementary 
path in Q with <71 = qn or qi,qn £ Vu and let s be its signature. For i £ |1,n — 1|, 
let ji be the index corresponding to the edge We have 


N 


^ S{j) / {f 23 -l{x) + f2j{x)) dx 

2 = 1 

N .Lj N 

= 2^s(j) / u){x)dx = 2'^s{j){uj{Lj)-Uj{Q)) 

n—1 

= ‘^^{up{q^+i) - Uj.iqi)) = 2 - Mji(9i)) = 0, 


i=l 


and thus / £ Yp{R). 

Conversely, take / £ Yp{R). For j £ |l,iV|, define Vj : [0,Lj] —>• C by 
„ -x)- f23ix) 

— o 


(41) 


One clearly has vj £ LP{[0,Lj],C). We define Uj as follows: let e £ f be the edge 
corresponding to the index j. Let (qi,..., qn) be any elementary path with qi £ Vu 
and qn = a{e). Let s : £ ^ {—1,0,1} be the signature of that path and, for 
i £ |1, n — 1], let ji be the index associated with the edge {< 7 ^, qi+i}. For x £ [0, Lj], 
set 


n—1 

lix) = '^s 

2=1 


(P) f 

Jo 


/2ji-l(C) + f2jM) 




f2j-l{Lj -£,) + f2j{0 


df. (42) 


This definition does not depend on the choice of the path [qi,... ,qn) with qi £ 
Vu and qn — Q:{e) thanks to the definition of the matrix R. It is an immediate 
consequence of (42) that {u,v) £ X^. The map / i-A {u,v) defines an operator 
S : Yp{R) —>• Xp. We clearly have T o S = IdYp(R) and S oT = Idx^, and thus T 
is bijective. The continuity of T and S follows immediately from (40), (41), and 
(42). □ 


Remark 5.4. When p = 2, one easily checks that -^T : X 2 —>■ Y 2 {R) is unitary. 

Remark 5.5. The operator T corresponds to the d’Alembert decomposition of the 
solutions of the one-dimensional wave equation into a pair of traveling waves moving 
in opposite directions. For every j £ |1, N}, f 2 j-i and f 2 j correspond to the waves 
moving from uj{j) to a{j) and from a{j) to uj{j), respectively (see Figure 1). 
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Figure 1. D’Alembert decomposition of the wave equation on the 
edge j G |l,iV]. 

Let us consider the operator B{t) in Yp(i?) defined by conjugation as 

Dm)) = {/ G Yp(A) I T-V G D(A(t))}, B{t)f = TA{i)mf. 

In order to give a more explicit formula for -B(t), we introduce the following nota¬ 
tions. 

Definition 5.6 (Inward and outward decompositions). The inward and outward 
decompositions of are defined respectively as the direct sums 

<?ev qev 

where, for every q G V, we set 

Wil = Span {{e 2 j \ w(j) = g} U {e 2 j-i \ a{j) = g}), 

W^out = Span {{e 2 j \ a{j) = g} U {e 2 j-i \ uj{j) = g}). 

For every g G V, we denote by 11?^ and the canonical projections of onto 
and respectively, which we identify with matrices in A4ng,2N(C). 

For n G N, let J„ denote the n x n matrix with all elements equal to 1. Set 
D = diag((—l)'^)jg|i_ 2 A]- For g G V and t G M, we set 

' (id„, n?„, if g G Vint, 

M«(t) = < ifgGVn, 

. i + vl{t) ^ 

We define the time-dependent matrix M = (wy )i,jG|i, 2 A] by 


M = -D 



D. 


(43) 


Remark 5.7. If the components of rj are nonnegative measurable functions, then 
M is measurable and its components take values in [—1,1]. 


Remark 5.8. Notice that and IF® are orthogonal whenever gi ^ g 2 , and 
similarly for the outward decomposition. Moreover, for each g G V, the spaces 
and are invariant under D. We finally notice that the image of M‘^{t) is 
contained in From these observations, we deduce that, for every g G V and 

t G R, 


mOMit) = -IllntMHt)D. 
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We finally obtain the following expression for B{t). 

Proposition 5.9. For t gM. and p G [1, +oo], the operator B{t) is given by 

{ 2N 

/£ Yp(i?)nnw^'’^([o,ir],c) 

27V 

m =Y.m,,{t)f,{L]), V* G [l,2iVl 

B{t)f = (45) 

Proof. Let / G Yp(i?) and {u,v) = T~^f G and notice that 

u'/x) = , v^{x) = ( 46 ) 

It follows from (40) and (46) that fi G W^’^([0, L)"], C) for every i G |1,2A^] if and 
only if Ui G W^’^([0, Li], C) and Vi G W^’P([0, Li], C) for every i G 

We suppose from now on that fi G W^’^([0, L)"], C) for every i G |1,27V|. Let 
Fq = (/*(0))*Gli,2Ari and Fl = {MLJ))zeli, 2 N}- The condition 

2N 

m = Y. m,,(t)/,(LJ), V* G [1,21V1 (47) 

i=i 

can be written as LJ) = M(t)FL. Thanks to the outward decomposition of 
this is equivalent to nQ^^DFo = Bl^^DM{t)FL for every q gV. By Remark 5.8, we 
have Il1^^DM(t) = —Ill^^M‘>{t)D, and thus (47) is equivalent to 

Ul^.DFo + Ul,,M'^it)DFL = 0, Vg G V. (48) 

If g G Vd, let j be the index corresponding to the unique edge in £q. To simplify 
the notations, we consider here the case a{j) = q, the other case being analogous. 
Then 

n«„,CFo + n‘‘„,M’{t)DF^ = 

=f2,m - = .'( 0 ) - „,( 0 ) - (.'( 0 ) + „,( 0 )) 

which shows that the left-hand side is equal to zero if and only if one has Vj{q) = 
—? 7 q(t)^(g). If g G Vu, the same argument shows that the left-hand side is equal 
to zero if and only if Vj{q) = 0. 

Finally, if q G Vint, one easily obtains that 

bIdfl = (g(g) - , n^nti^Fo = (-$(g) - • 
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Since one has 

Ul^.DFo + nl,,M\t)DFL 

= (^-2vj{q) - ^ J2keEg (^( 5 ) “ 


I^Tx Jn 

Tin 


(S(9)-^4(9)) 


iG£g 


iG£g 


The right-hand side is equal to zero if and only if Vj{q) = Vk{q) for every j, k G £q 
and Efegf, ^( 9 ) = 0 . 

Collecting all the equivalences corresponding to the identities in (48), we conclude 
that (44) holds. 

Let now / G D[B{t)) and denote {u,v) = T~^f G D{A{t)), g = Then 

g = TA{t)T-^f = TAit)iu, v) = T{v, u"), 
and so, by (40), for every j G |1, 2iV|, 

52t-i(a;) = v'jiLj - x) + u"{Lj - x) 
d 
dx 

92j{x) = v'^(x) - u”{x) = {vj{x) - u'^{x)) = -f2j{x) 
which shows that (45) holds. □ 

The operator T : —>■ Yp(i?) transforms (39) into 

F{t) = B{t)F{t). 

This evolution equation corresponds to the system of transport equations 


= -x)+ u'j{Lj - x)) = -f2j-iix), 

A. 

dx 




i G |1, 2A^], t G [0, -boo), X G [0, LI], 

2N (49) 

= iG [l,2Af], tG [0,-boo), 

4=1 

where F{t) = (/i(t))ig|i, 2 A]■ The following property of the matrix M{t) will be 
useful in the sequel. 

Lemma 5.10. For every t G K, 

M(t)TM(t) = id2A- E 

Proof. Notice that, for every g G V, can be written as 

M^t) = (n^jT (^A,(t)id„,n^, 

where Xq{t) = if g G Vd and Xq{t) = 1 otherwise, while d, = 1 if g G Vint and 

Sq = 0 otherwise. By a straightforward computation, one verifies that, for every 

ge V, 

T 


Ag(t)ldnq ^ ^qJuq 


Xq{t)ldn^ - SqJn^ = Xq{t)'^ Un^ . 
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Noticing furthermore that, for every gi ,(;2 G V, one 

deduces that 


= D 


(JGV 


D. 


Since the term between brackets in the above equation is diagonal and Xq{t)^ = 
1 — for q G Vd, the conclusion follows. □ 


5.2. Existence of solutions. Thanks to the operator T : — >■ Yp{R), one can 

give the following definition for solutions of (38). 


Definition 5.11. Let Uo G Xp and r] = {riq)qeVi be a measurable function with 
nonnegative components. We say that U : M+ —>■ X^ is a solution of T,i^{0,L,r]) 
with initial condition Uq if T~^U : R+ —^ Yp(i?) is a solution of (49) with initial 
condition T~^Uo G Yp(i?). 

For every Fq G Yp(i?), it follows from Proposition 4.3 that (49) admits a unique 
solution F : IR_|_ —)> Xp. In order to show that this solution remains in Yp(ii) for 
every t > 0, one needs to show that Yp(i?) is invariant under the flow of (49). 


Proposition 5.12. For every t G R, RM{t) = R. 

Proof. Thanks to the inward decomposition of we prove the proposition by 
showing that for every q G V and t G M, 


- RD{Ul^,y 


Aq(t)Idng __ SqJrig 


= Rmi 


(50) 


where Xq{t) and 6q are defined as in the proof of Lemma 5.10. Without loss of 
generality, it is enough to consider the case where i? is a line matrix, i.e., we consider 
a single elementary path {qi,... ,qn) in t/ with qi = qn or qi,qn G Vu, with signature 
s. Then R = (pj)j 6 |i, 2 Ari is given by p 2 j-i = P 2 j = s{j) for j G |l,iV|. For 
i G |l,n — 1], denote by ji the edge corresponding to {qi,qi+i}. Let us write 
R = + e 2 ji)"'" and notice that 

n—1 

RD = Y^ s{ji){-e2j,-i + e2jj^. 

By definition of the signature s, one has, for i € |1, n — 1|, 

- (nS)^n-], 

= e?,. [(nrrnr-(n-)^n-], 

and 

[(n^;t)^n^;t - (nrt^)^nrtl • 

One deduces that 

n—1 

RD = E(e2,.-1 + e2,,r - (n:(;)^ni(;] 


n—1 


= E(e2..-1 + 62 ,.)^ [(n^‘ut)^n^;t - (nrt^)^nrt^] • 
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By using the above relations, Equation (50) can be rewritten as 


n—1 

^out ~ + ^2ji) 

n—1 

= (^OTi+i ~ ^qqi) + ^2ji)- (51) 

i=l 

Such an identity is trivially satisfied if q ^ {gi,..., q„}. Assume now that either 
q = qi for some i € | 2 ,u — 1 ] or g = gi = g„ (and in the latter case set i = n and 
define jn+i = ji)- In particular, q € Vint and Xq{t) = 5q = 1. We therefore must 
prove that 


Ag(t)Idng ^qJriq 

Tlq 




1 + e2ji_i - e2ji-i - e2ji) 

= II?n(e 2 ii_i-i + e 2 ji_i - e2ji-i - e2ji)- (52) 


By definition of B?) and , one has that 

nout(e 2 ii_i-i + e 2 ji_i - e 2 j,-i - e 2 jJ = n?;(e 2 j,_i-i + e 2 ji_^ - e 2 ji-i - e 2 ji) = w, 


where w € has all its coordinates equal to zero, except two of them, one equal 
to 1 and the other one equal to —1. Hence Jn^.w = 0 and (52) holds true. 

It remains to treat the case g S {gi, g„} C Vu- In this case, Xq{t) = 1 and 6q = 0, 
and we furthermore assume, with no loss of generality, that g = gi. We can rewrite 
(51) as 

nout(e2ji-i + e2ji) = n?((e2ji-i + 62^1), 

which holds true by definition of H?) and This concludes the proof of the 

proposition. □ 


The main result of the section, given next, follows immediately from Propositions 
4.5 and 5.12. 

Proposition 5.13. Let {G,L) be a network, p € [l,+oo], and ij = (r]q)q^y^ be a 
measurable function with nonnegative components. Then, for every Uq G Xf, the 
system Tii^(G,L,ri) defined in (38) admits a unique solution U : M+^X-. 

5.3. Stability of solutions. We next provide an appropriate definition of expo¬ 
nential type for (38). 

Definition 5.14. Let {Q, L) be a network, p G [1, -boo], and H be a subset of the 
space of measurable functions 77 = (qq)qeVd with nonnegative components. Denote 
by Sdj(G, L, V) the family of systems E^(t/, L, ij) for ij gV. We say that S,^(t/, L, V) 
is of exponential type 7 in if, for every e > 0, there exists K > 0 such that, for 
every ij G T) and uq G X^, the corresponding solution u of Tiu:{G, L, q) satisfies, for 
every t > 0, 

||R(t)||p<i^e(^+^)‘|K||^. 

We say that L, D) is exponentially stable in X^ if it is of negative exponential 

type. 

Given D as in the above definition, we define 

M. = {M : R —>• Al 2 Ar(R) I XL is given by (43) for some 77 G V}. 
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Thanks to the continuity of T and T~^ established in Proposition 5.3, we remark 
that is of exponential type 7 in if and only if is of expo¬ 

nential type 7 in Yp(i?). As a consequence of Corollary 4.11, we have the following 
result in the case of arbitrarily switching dampings rjq, q € Vd- 

Corollary 5.15. Let (1/,A) be a network, d = #Vd, S a subset of and 

D = The following statements are equivalent. 

i. Yii^{Q,K,'D) is exponentially stable in Xf for some p G [l,-|-oo]. 
ii. L, V) is exponentially stable in for every L G V+(A) andp G [1, -foo]. 

We can now provide a necessary and sufficient condition on and S) for the 
exponential stability of A, L°°(R, S))). 

Theorem 5.16. Let (f/, A) be a network, d = #Vd, 2) a bounded subset of 
and T) — S)). Then E,^(5,A,21) is exponentially stable in X^ for some 

p G [1, -foo] if and only if Q is a tree, Vu contains only one vertex, and D C 

Proof. Similarly to Remark 4.12, the exponential stability of A, V) is equiva¬ 
lent to that of A, 23)). We therefore assume with no loss of generality 

that 23 is compact. 

Suppose that either G is not a tree, Vu contains more than one vertex, or 23 
contains a point rj with rj.^ — 0 for some g S Vd. Let (gi,..., g„) be an elementary 
path in G with gi = gn > gi > g™ G Vu, or gi G Vu and qn = q. Let s be its signature 
and, for i G |1,u — 1|, let ji be the index corresponding to the edge {g^, gi+i}. Take 
L G V+(A) nN^, which is possible thanks to Proposition 3.9. For j G |1, A], we 
define 



s{ji) sin(27rt) sin(27ra;), if j = ji for a certain i G |1, n — 1], 
0 , otherwise. 


One easily checks that (wj)jG[i,Ail is a solution of Tii^{G, L,r]) for every p G T). 
Since it is periodic and nonzero, L,23) is not exponentially stable in X^ for 

any p G [1, -foo], and so, by Corollary 5.15, T^ujiG, A, 23) is not exponentially stable 
in Xp for any p G[1, -foo]. 

Suppose now that t/ is a tree, Vu contains only one vertex, and D = D Gl (R;!^)'^. 
Up to changing the orientation of G, we assume that a(j) = g for every g S Vu 
and j G £q. Let 77nun = rnin^gj) minggVd’ll? > 0- Let U = {u,v) be a solution 


of A,23) in X 2 and / = -^TU. Notice that ||/( 0 IIy 2 (R) ~ I|t^(^)ll 2 thanks 

to Remark 5.4. For t > 0, denote Fb(t) = (/i(L 0))ig|i 2 Ar] and Fa(0) = {fi{t, 
2Ar|, SO that Fb(t) = M(t)F\{f). For t > 0 and s G [0,Amin], we have, by 
Lemma 5.10, 
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= iit'(f)ii2 -1: i: 


rt-\-s 


4?7g(T) 


gGVd i&£q 


't {^ + Vqir)y 


|/ 2 i-i(r, Ai)^ dr, 


and, since 

/2j-i(t, Ai) = 
we conclude that 


+Vi{T,0) _ l + ? 7 ,(r) duq 


\/2 dx 

\\u{t + s)\\l = \\uml- E E r'^Vqir) 

9GVd ie£q 


(r,0), Vg e Vd, Vi e f,, 




dr. (53) 


Since (53) holds for every t > 0 and every s G [0, Amin], one can easily obtain by an 
inductive argument that it holds for every t > 0 and s > 0. Hence, for every t > 0 
and s > 0, 


\mt + s)\\i 



dui 

dx 


("^, 0 ) 


2 

dr. 


We can thus proceed as in [11, Chapter 4, Section 4.1] (see also [26]) to obtain the 
following observability inequality: there exist c > 0 and £ > 0 such that, for every 
t > 0, 





2 

dT>c||;7(< + £)|[^ 


This yields the desired exponential convergence in X 2 , and hence in for every 
p G [1, -\-oo] thanks to Corollary 5.15. □ 
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